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We use the high-order nonlinear Schrodinger equation (NLSE) derived to model the evolution 
of slowly modulated wave trains with narrow spectrum on the surface of ideal finite-depth fluid. 

This equation is the finite-depth counterpart of celebrated Dystlie’s equation, which is usually used 
for the same purpose in the case of infinite depth. We demonstrate that this generalized equation 
admits bright soliton solutions for depths below the modulation instability threshold kh « 1.363 (fc 
being the carrier wave number and h the undisturbed fluid depth), which is not possible in the case 
of standard NLSE. These bright solitons can exist along with the dark solitons that have recently 
been observed in a water wave tank [Phys. Rev. Lett. 110, 124101 (2013)]. 

PACS numbers: 05.45.Yv, 47.35.Bb, 92.40.Qk 


Bright and dark solitons are the fundamental self- 
localized modes of the optical field in nonlinear disper¬ 
sive media such as waveguides, optical fibers, and pho¬ 
tonic crystals [l] . Bright solitons are characterized by a 
localized intensity peak on a homogeneous background, 
while dark solitons can be described by a localized inten¬ 
sity hole on a continuous wave background [2]. Solitons 
are known to exist due to the balance of dispersion and 
nonlinearity and propagate without changing their shape 
and keeping their energy [3]. Bright solitons are formed 
when the group-velocity dispersion in an optical fiber is 
anomalous (or, similarly, when the nonlinearity of a pla¬ 
nar waveguide is self-focusing). In this case, the uniform 
carrier wave is unstable with respect to long-wave mod¬ 
ulations allowing for the formation of solitons. This type 
of instability is known as the modulation instability [4|. 
On the contrary, dark solitons are formed in the case 
of normal group-velocity dispersion in fibers (or a self- 
defocusing nonlinearity in waveguides), when a uniform 
carrier wave is modulationally stable. 

In the case of water waves, bright solitons are known 
to appear in the form of surface envelopes of modulated 
wave trains when the uniform carrier wave is modula¬ 
tionally unstable [sj. This happens for water depths h 
above the modulation instability threshold, namely, at 
kh > 1.363, k being the carrier wave number. In addi¬ 
tion to theoretical predictions, envelope solitons were ob¬ 
served experimentally in Refs. , mostly in the case of 
deep water (kh 1). Dark solitons can appear on shal¬ 
low water below the modulation instability threshold, at 
kh < 1.363 0 - They have recently been observed in 
a series of experiments performed in a water wave tank 

mm- 

In mathematical terms, bright and dark solitons are de¬ 
scribed by the nonlinear Schrodinger equation (NLSE) of 
the focusing and defocusing types, respectively 0S0. 
NLSE takes into account the second-order dispersion and 
the phase self-modulation (cubic nonlinear term). In the 
general context of weakly nonlinear dispersive waves, this 


equation was first discussed by Benney and Newell la¬ 
in the case of gravity waves propagating on the surface 
of infinite-depth irrotational, inviscid, and incompress¬ 
ible fluid, NLSE was first derived by Zakharov [lij]. The 
finite-depth NLSE was first derived by Hasimoto and 
Ono 0. 

To achieve a better comparison with experiment and 
to trace nonlinear effects in the vicinity of the modula¬ 
tion instability threshold kh ~ 1.363, where the lead¬ 
ing nonlinear term vanishes, high-order nonlinear and 
nonlinear-dispersive effects should be taken into consid¬ 
eration. In the case of infinite depth, such a hi gh-o rder 
NLSE (HONLSE) was first derived by Dysthe 0 and 
then rewritten in a more general form by Trulsen et al. 
0 . It includes the high-order dispersion and cubic non¬ 
linear dispersive terms as well as an additional nonlinear 
dispersive term describing the input of the wave-induced 
mean flow. Dysthe’s equation was extensively used in nu¬ 
merical simulations of wave evolution and showed good 
agreement with experiment and_ simulations based on 
fully nonlinear Euler equations L9j-24|. 

In the case of finite depth, HONLSE is generally cou¬ 
pled with additional equation for the wave-induced mean 


flow [251,[26|. No modeling of water wave evolution has 
been performed with these general equations because of 
their complexity. Sedletsky [27] used an additional power 
expansion of the induced mean flow to derive a single 
HONLSE for the first-harmonic envelope of surface pro¬ 
file. This equation is the direct counterpart of Dysthe’s 
equation 17, 18] but for the case of finite depth. Slun- 
yaev 28] confirmed the results obtained in Ref. 27| and 
extended them to the next order. Gandzha et al. [29] 
rewrote this equation in dimensionless form and used it 
to model the evolution of bright solitons on finite depths. 
When the HONLSE terms were taken into considera¬ 
tion, one-soliton solutions to NLSE were transformed 
into quasi-soliton solutions with slowly varying ampli¬ 
tude. These quasi-solitons were found to propagate with 
nearly constant speed and possess the unique property of 
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solitons to exist over long periods of time without break¬ 
ing. Their speed was found to be higher than the speed 
of the bright NLSE solitons taken as initial conditions 
in computations. This phenomenon was observed earlier 
both in experiment and numerical modeling in the case 
of deep-water limit in Refs. 20j, 30]. 

In the present work we demonstrate that the HONLSE 
derived in Ref. [HJ admits exact solutions in the form of 
bright solitons below the modulation instability thresh¬ 
old. Such bright solitons can be observed in the same ex¬ 
perimental setup that was used by Chabchoub et al. 11] 
to observe dark solitons at kh = 1 . 2 . 

We start from the fully nonlinear Euler equations writ¬ 
ten for the potential two-dimensional waves on the sur¬ 
face of irrotational, inviscid, and incompressible fluid un¬ 
der the influence of gravity. Waves are assumed to prop¬ 
agate along the horizontal 2 -axis (—oo < x < oo), and 
the direction of the vertical y-axis is selected opposite to 
the gravity force. The fluid is assumed to be bounded 
by a solid flat bed y = — h at the bottom and a free sur¬ 
face y = 77 ( 2 , t ) at the top. The atmospheric pressure 
is assumed to be constant on the free surface. Then the 
evolution of waves and associated fluid flows is governed 
by the following set of equations (3l1 ]: 


*&xx ^yy — O 7 

< 77 ( 2 , t); 

( 1 ) 

+ + +gy = 0, 

y = v(x, t); 

( 2 ) 

r)t-$y + rjx^x = 0 , 

y = v(x, ty, 

(3) 

II 

O 

-d 

1 

II 
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Here, < 1 >( 2 , y , t ) is the velocity potential, g is accelera¬ 
tion due to gravity, t is time, and the indices designate 
the partial derivatives over the corresponding variables. 
The linear dispersion relation for wave trains with carrier 
frequency w and wave number k is 

uj 2 = gktanh(kh) = gka, a = tanh (kh), (5) 


and the carrier group speed is 

v ’ = W; = lLy-'’ 2)kh + °)- <6) 

Considering slowly modulated wave trains with narrow 
spectrum around the carrier frequency and wave num¬ 
ber, the unknown free-surface displacement and velocity 
potential can be looked for in the form of Fourier series 
with slowly variable coefficients: 



exp(in(fca: 
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* 
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- vt)), 
(7) 


where * stands for complex conjugate. In the small- 
amplitude approximation the unknown functions r) n and 
can be all expressed in terms of the first-harmonic 


amplitude iji with the use of additional multi-scale power 
expansions 27]. In this way, the original set of equations 


<[T|) (HI) can be reduced to one evolution equation for the 
complex-valued amplitude 771(2 , t). Introducing dimen¬ 
sionless time, coordinate, and amplitude 


w 

r=-t , x = kx, 
c 


= 2y/ckrji, 


( 8 ) 


this evolution equation can be written in the following 
form [29i ]: 


u T = —aiu x — ici2U xx + iao : o, 0 \u\ 2 u 

+ ( a 3»xxx - ai.o.o^xM 2 _ ao,o,i«X)’ ( 9 ) 

This HONLSE describes the evolution of the first- 
harmonic envelope of the surface profile with taking into 
account the third-order dispersion and cubic nonlinear 
dispersive terms. Here, c and ai are the dimensionless 
phase and group speeds, respectively, 

4 a 2 

c= ->0, (10) 

v 

01 = — Vg = — — ^(l — a 2 ) kh + o'j >0, (11) 

v = (a 2 — l) (3cr 2 + l) k 2 h 2 — 2a (a 2 — l) kh — a 2 . 
The second- and third-order dispersion coefficients are 


“ 2 - 2 ’ 
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,( 1 ) 


a 
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= -3a 2 (a 2 - 1) 
= —3a (a 2 — l) (3a 2 + l) , 

= (a 2 - 1) (15a 4 - 2a 2 + 3) , 


and the coefficient at the cubic nonlinear term is 

ao,o,o = -Y&p^Yj a oXo( kh ) P ’ 


p =0 

.4 


a. 


4, 0,0 = ~° 2 (7o 4 - 38a 2 - 9) , 

,£4 0 = 2a (3a 6 - 23a 4 + 13a 2 - 9) , 

a o%,o = (° 2 - 4 2 ( 9(j4 - 10<j2 + 9 ) > 
v = (a 2 — l) 2 k 2 h 2 — 2a (a 2 + l) kh + a 

The cubic nonlinear dispersion coefficients are 
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0 = a 5 (a 6 - 40a 4 + 193a 2 + 54) , 


a, 


(0) 0)1 = -a 5 (a 6 - 7a 4 + 7a 2 - 9) , 


4(4 0 = -a 4 (a 8 - 109a 6 + 517a 4 + 217a 2 + 270) 
4(o, 1 = o 4 (a 8 + 20a 6 - 158a 4 - 28a 2 - 27) , 


( 12 ) 


(13) 
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FIG. 1. HONLSE coefficients as functions of kh. 


2 (l — a 2 ) kh + a 1 3 — a 2 

&o = -> oi i = a.2 = — —o-, 

cv tJc oca” 3 

where 6 = kx — cot = x ~ CT is the wave phase. The 
similar_expression for the velocity potential is given in 
Ref. [H. 

The region of modulation instability of a homogeneous 
solution u(x, t) = uq exp (iao, o, o|'“o| 2 ' r ) to Eq. © can 
be determined by introducing a small perturbation to the 
amplitude uq (as described, e.g., in Refs. 

«(x, T ) = (“o + e(x, r)) exp (i“o,o,oK| 2 T) , (16) 

e(x> T ) = 4 ex P ( iK X - + eg exp (-irex + ifl*r ). 


a^o 0 = —2a- 3 (3a- 10 + 18a- 8 — 146a- 6 
-172a- 4 + 183a- 2 - 270), 
a[, 2) 0i 1 = 2a- 3 (3a- 10 - 63a- 8 + 314a- 6 
—218a- 4 + 19a- 2 + 9), 

a i 3 o,o = 2a 2 (“ 2 - t) (7a- 10 - 58a 8 + 38a 6 
+52a 4 -181a 2 + 270), 
a ( 0 3 { 1 = -2a 2 (a 2 - 1) (7a 10 - 79a 8 + 282a 6 
-154a 4 -a 2 + 9), 

a ito,o = a (1 - a 2 ) 3 (lla 8 - 99a 6 - 61a 4 + 7a 2 + 270) , 
o,q 4 q ! — a (a 2 — l) 3 (lla 8 — 48a 6 + 66a 4 + 8a 2 + 27) , 
a po,o = (o' 2 - l) 5 (3a 6 - 20a 4 - 21a 2 + 54) , 
a$,l = — (o' 2 — l) 5 (3a 6 + 7a 4 — lla 2 + 9) . 


The coefficients 03, a0,0,o, “1,0,0, and “o, 0,1 are all real 
and depend only on one dimensionless parameter kh. 
Their behavior as functions of kh is shown in Fig. |T| 
It can be seen that Eq. © is valid at kh > 1 , where its 
coefficients do not diverge. At kh —> 00, the following 
asymptotics are easily obtained: 

1 13 1 

“3 — —, ao,o,o — — 2’ a bo,o — 2’ “0,0,1 —-p 


They coincide with the corresponding coefficients of Dys- 
the’s equation 18 ], except for the term including the 
wave-induced mean flow. This term cannot be explic¬ 
itly reconstructed from Eq. © because of the addi¬ 
tional power expansion of the wave-induced mean flow 
made to derive this equation. However, it can be re¬ 
constructed from the equations generating Eq. © , at 
the stage when the wave-induced mean flow has not 
been excluded from the equation for u [ 27 |. By setting 
“3 = “1,0,0 = “0,0,1 = 0 , Eq. © is reduced to the stan¬ 
dard NLSE derived by Hasimoto and Ono 0 - 

The dimensionless free surface displacement is ex¬ 
pressed in terms of u as follows 


kr) = otu\u \ 2 + at\ Re(uexp(i 0 )) + 2 ci 9 Re(u 2 exp( 2 i 0 )), 

( 15 ) 


Here, we assume the perturbation frequency Q to be 
complex-valued and the perturbation wave number n to 
be real. Substituting this ansatz in Eq. © leads to the 
following dispersion relation between Q and re: 

ft = (ai + «i,o,o|“o| 2 ) « + a 3 n 3 

± reY /, 2a 2 “o,o,o|“o| 2 + “o,o,il w o| 4 + “|k 2 - (17) 

A homogeneous solution is modulationally unstable when 
the perturbation exponentially grows with time. This 
happens when Im ft > 0 , which effectively requires the 
radicand in Eq. G 3 to be negative. This condition is 
satisfied when 


«o,o,o <—“0,0,il“o| 2 , ( 18 ) 

where we took into account that “2 = ^- In the NLSE 
case, when 00,0,1 = 0, condition m is reduced to the 
well-known modulation instability criterion ao,o,o < 0, 
which holds true at kh > 1.363. In the HONLSE case, 
this threshold slightly shifts to higher kh , depending on 
the value of uq. 

Bright and dark solitons are defined as [3] 

ub(Xi t) = u 0 sech(A'(x - W)) exp (irex - ifI t) , (19) 
ud(x, t) = uq tanh(Ar(x — Hr)) exp (irex — iflr). (20) 

Here, u 0, re, ft, and V, are the soliton’s complex ampli¬ 
tude, wave number, frequency, and speed, respectively. 
The soliton’s initial position was chosen to be located at 
X = 0. In the NLSE case, bright one-soliton solutions 
of form m were first derived in Ref. 32]. The follow¬ 
ing relationships between the soliton parameters can be 
established in this case: 

K = |uo|x/-“o,o,0) 

ft = reai + ^ (K 2 ~ ft 2 ) , V = ai — re. (21) 

Bright NLSE solitons can exist at a 0,0,0 < 0, which ex¬ 
actly corresponds to the region where a homogeneous so¬ 
lution is modulationally unstable. The parameters uq 
and re are free parameters of the problem. However, of 
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FIG. 2. Wave number re of the HONLSE soliton versus kh. 


physical relevance are only such amplitudes and wave 
numbers that preserve the conditions of small amplitude 
and narrow spectrum: |uq| -C 1 , re <C 1. 

Dark solitons of form (l20l) are referred to as fundamen¬ 
tal dark solitons 0. In the NLSE context^ they were first 
derived in a more general form in Ref. |33|. Dark NLSE 
solitons exist below the modulation instability threshold, 
at ao,o,o > 0, and the corresponding relationships be¬ 
tween the parameters are 


I \ — |uq|\/oo,o,Oj 


fi = rear - - (2 K 2 


+ re 2 ), V — ai — 


re. 


( 22 ) 


Relations m and ( 1221 ) are no longer valid in the 
HONLSE case, when the high-order terms are not set 
equal to zero. In this case, bright NLSE solitons are 
transformed into quasi-solitons, which were obtained by 
numerical integration of Eq. (©) in Ref. (29). Below we 
prove that HONLSE (©) admits a new family of exact 
bright soliton solutions that exist below the modulation 
instability threshold kh ~ 1.363. Substituting ansatz 
(fTiJl) in Eq. ©, one can obtain the following relations: 


K = \u 0 \VS, s = - ai ’°’° r +ao ' °' 1 


re = 


6a 3 

ai,o,o + °o, o, l — 60300,0,0 


(23a) 

(23b) 


n = 


12a 3 ao,o, 1 

rear + ^ 2 (1 — 6rea 3 ) — re 2 (1 — 2rea 3 )^, (23c) 

V = ai - re + (3re 2 - K 2 ) o 3 . (23d) 


Note that such relations were derived earlier for the gen¬ 
eral form of Eq. © in Refs. 34 -HI without referring to 
water waves. Anzats csd is an exact solution to Eq. m 
when the radicand in Eq. (I23a[) is positive, S > 0. This 
condition holds true in the following depth range: 



FIG. 3. Bright and dark solitons of HONLSE © for kh = 1.2 
and uo = 0.17. The bright soliton parameters are determined 
by relations (1231) : K « 0.052, re « 0.23, fl f» 0.30, V 
1.18. The dark soliton parameters are chosen such that they 
closely correspond to the water-tank experiments described 
in Ref. Q: Kk 0.066, re = 0, fi « -0.0043, V = ai « 1.39. 


which lies below the modulation instability threshold. 
The soliton wave number re is not a free parameter any 
longer. Its behavior as a function of kh is shown in Fig. [2] 
The condition of narrow spectrum (slow modulation) 
holds true only in a narrow range around kh « 1.249, 
where re = 0. At kh = 1.2, we have re « 0.23. This wave 
number can still be regarded small as compared to unity, 
but not so small as in the NLSE case. Figure [3] shows 
a bright soliton solution to HONLSE ©) along with the 
dark soliton that exists at the same dimensionless depth 
kh = 1.2. The corresponding free surface profile com¬ 
puted by formula (1151) is shown in Fig. [I] The parame¬ 
ters of the dark soliton shown in Fig.[3]were selected such 
that they closely correspond to the water-tank geome¬ 
try and dark-soliton experiments described in Ref. Q 
These experiments were performed at the same depth 
(kh = 1.2), the carrier wave number was k = 3 m — , 
and the carrier amplitude was a = 0.04 m. The corre¬ 
sponding dimensionless amplitude ka = 0.12 corresponds 
to the amplitude of soliton shown in Fig. 01 The evo¬ 
lution time t = 27.6 (t = 10.8 s) was chosen such that 
the wave group travels the distance L = V g t = 12.8 nr 
(V g = 1.19 m/s) to the last wave gauge position along 


the wave tank in the experimental setup 11 . 


0.763 < kh < 1.222, 


(24) 


The dark soliton shown in Fig. [3] was computed by in¬ 
tegrating HONLSE (©) numerically using the split-step 
Fourier technique described in [29J with the dark NLSE 
soliton taken as the initial wave form. On the adopted 
time scales, the corresponding numerical HONLSE so¬ 
lution identically followed the initial dark soliton shape. 
Deviations from NLSE can be observed on much larger 
time scales (r > 2000) in the form of slight oscillations 
of the wave envelope. Besides such dark quasi-solitons, 
HONLSE (©) admits a family of exact dark soliton solu- 
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FIG. 4. Free surface profile with envelope corresponding to 
the bright soliton shown in Fig. \3\(kh = 1.2). 


tions defined by the following set of parameters: 


K=\u 0 \V^S, (25a) 

k 2 

= Kcii + K 2 (6ftci3 — 1)-— (1 — 2na 3) , (25b) 

V = a± — k + (3k 2 + 2A' 2 ) a 3 , (25c) 


with k given by the same formula (123bl) as for the bright 
solitons. These relations can be obtained in the same way 
as relations (l23l) . Similar expressions for dark HONLSE 
solitons were also derived in Ref. (34| without referring 
to water waves. Dark HONLSE solitons exist at S < 0, 
i.e., at kh > 1.222. In contrast to bright solitons, NLSE 
and HONLSE dark solitons can exist in the overlapping 
depth ranges. 

Thus, we presented for the first time new families of 
bright and dark solitons that exist below the modula¬ 
tion instability threshold for water waves on the surface 
of ideal finite-depth fluid. These solitons are the exact 
solutions to the HONLSE derived earlier in Ref. [ 2 ^] to 
describe the evolution of slowly modulated wave trains. 
The results of this study can readily be used in water- 
tank experiments to supplement the recent observations 
of dark solitons described in Ref. [11]. 

We would like to thank Dr. Denys Dutykh for helpful 
discussions. 
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